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ABSTRACT 

Dust polarization observational results are analyzed for the high-mass star for- 
mation region W51 from the largest parent cloud (~ 2 pc, JCMT) to the large-scale 
envelope (~ 0.5 pc, BIMA) down to the collapsing core e2 (~ 60 mpc, SMA). Mag- 
netic field and dust emission gradient orientations reveal a correlation which becomes 
increasingly more tight with hig her resolution . The previously developed polarization 



intensity gradient method (|Koch et al.ll2012|) is applied in order to quantify the mag 



netic field significance. This technique provides a way to estimate the local magnetic 
field force compared to gravity without the need of any mass or field strength measure- 
ments, solely making use of measured angles which reflect the geometrical imprint of 
the various forces. All three data sets clearly show regions with distinct features in 
the field-to-gravity force ratio. Azimuthally averaged radial profiles of this force ratio 
reveal a transition from a field dominance at larger distances to a gravity dominance 
closer to the emission peaks. Normalizing these profiles to a characteristic core scale 
points toward self- similarity. Furthermore, the polarization intensity-gradient method 
is linked to the mass-to-flux ratio, providing a new approach to estimate the latter one 
without mass and field strength inputs. A transition from a magnetically supercritical 
to a subcritical state as a function of distance from the emission peak is found for the e2 
core. Finally, based on the measured radius -dependent field-to-gravity force ratio we 
derive a modified star formation efficiency with a diluted gravity force. Compared to a 
standard (free-fall) efficiency, the observed field is capable of reducing the efficiency 
down to 10% or less. 



1 Academia Sinica, Institute of Astronomy and Astrophysics, Taipei, Taiwan 

2 Universite de Bordeaux, Observatoire Aquitain des Sciences de l'Univers, 2 rue de FObservatoire, BP 89, F-33271 
Floirac Cedex, France 

3 CNRS, UMR 5804, Laboratoire d'Astrophysique de Bordeaux, 2 rue de 1' Observatoire, BP 89, F-33271 Floirac 
Cedex, France 

4 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 



-2- 



Subject headings: ISM: clouds — ISM: magnetic fields, polarization — ISM: individ- 
ual (W51, W51 e2) — Methods: polarization 



Introduction 



The magnetic field is being recognized as a crucial component in the star formation process. 
Evidence for its presence and significance is growing with the advance of an increasing number 
of instruments capable of providing high-quality polarization observations. To date, various ob- 
servations - based on different emission or absorption mechanisms over a range of wavelengths 
- cover scales from a few pc down to mpc in molecular clouds. Nevertheless, the precise role of 
the magnetic field and its interplay with, e.g. turbulence and gravity, remain a debated topic in 
the current literature. In order to make further progress, knowledge of both the field morphology 
and the field strength needs to be combined. It is unfortunate that many of the currently avail- 
able observational techniques can only address one or the other of these field properties. Whereas 
dust polarization reveals a plane of sky projected magnetic field orientation and morphology (e.g., 
Hildebr and! 1 1 9 8 8|) . it does not provide any inform ation about the magnetic field strength. On the 
other hand, Zeeman splitting in spectral lines (e.g. JCrutcher et al.ll2009|) gives a line-of-sight field 
strength, but typically the field morpholo gy is not recovered due t o the only isolated detect ions. An 
exception might be the recent works by ISurcis et al.l (|201l|) and Iviemmings et al.l (|2011r) . where 
the morphology starts to become visible based on several tens of maser detections. 

Besides the intrinsic shortcomings of the observing techniques mentioned above, analysis 
tools and interpretation of polarization data are still in a developmental stage. Consequently, the 
benefit of the additional polarization information beyond its mere detection and imaging pur- 
pose is often not very obvious. With the growing number of high-quality polarization obser- 
vations from various instruments, it is, thus, paramount to investigate new methods and strate- 
gies to further explore and reveal the physics hidden in polarization data. In order to make fur- 
ther progress, efforts on several fronts are needed and eventually must be combined. On one 
hand, magneto-hydrodynamic simulations producing synthetic observati onal maps can charac 



terize generic features and provide guidelines for data interpretation (e .g., iNakamura & Lill201 1 



Falceta-Goncalves et al.ll2008l : lli & Nakamura! 12004 1 Allen et al.l 120031) . It is further desirable to 
include radiative transfer modeling in this approach. On the other hand, methods and techniques 
inspired by observed data can lead to new unexpected insights in a phenomenological way. Such 
recent approach es are, e.g., the polarization dispersio n function used to trace the relative turbu- 
lence level (|Houde et al.ll2009i : iHildebrand et al.ll2009|) a nd the ambipolar dif f usion scale isolate d 
by comparing coexisting ion and molecular line spectra (|Hezareh et al.ll201 It Li & Houdell2008|) . 
Here, we are further exploring the recently developed polarization - intensity gradient method 
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(IKoch et al.ll2012|) . This technique leads to a local position-dependent measurement of the mag- 
netic field strength, and it additionally provides an estimate of the local field-to-gravity force ratio 
in a model-independent way. 

We are applying our method to a set of dust emission polarization data in the sub-millimeter 
regime, where the dust grains are though t to be aligned with their sho r ter ax i s parallel to the mag- 



netic field lines due to radiative torques (lDraine & Weingartnerll 19961 1 19971 : lLazarian 



emitt ed light, therefore, appears to be polarized perpendicular to the field lines (e.g., 



2000). The 



Hildebrand 



1988|) . Dust polarization observations are now bein g routinely carried out w it h the Submillime 



Rao et a 



(2009); 



ter Array (SMA) (Tang e t al. 2012; in preparation; iTang et all (I2010[l2009bl a|): 
Girart et all fowA I2006I) ). This study is part of the program on the SMaR (Iho et al.1 l2Q04h to 
investigate the structure of the magnetic field from large to small scales. 

The paper is organized as follows. Focusing on the magnetic field detections in the W51 star 
formation region, Section[2]describes the relevant polarization observations fro m the JCM T, BIMA 
and the SMA. The key results of the polarization - intensity gradient method (IKoch et al.ll2012|) - 
which serve as a starting point for this work here - are summarized in Section [3] Section |4] starts 
with pointing out the correlation between the magnetic field and intensity gradient orientations, 
and then presents our results on the magnetic field significance. Implications of our findings for 
the mass-to-flux ratio and the star formation efficiency are discussed in Section[5l A summary and 
conclusion are given in Section [61 



2. Observations and Source Description 

The W51A| cloud is one of the most act ive and luminous h igh mass star formation sites in 
the Galaxy. Located at a d istance of ~7 kpc (|Genzel et al.N 198 lh . 1" is equivalent to ~ 0.03 pc. 
Chrysostomou et al.l (|2002|) measured the polarization at 850/im with SCUBA on the JCMT across 
the molecular cloud at a scale of 5 pc with a binned resolution 9 ~ 9."3. Their observation encom- 
passes two cores with polarized emission detected both in the cores and the region in between them. 
At this scale, the polarization appears to be organized but with a morphology that changes from 
the dense cores to the surrounding more diffuse areas. lLai et al.l (|200l|) reported a higher angular 
resolution (3") polarization map at 1.3 mm obtained with BIMA, which resolved out large-scale 
structures. In contrast to the larger scale JCMT map, the polarization appears to be uniform across 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics, and is funded by the Smithsonian Institution and the Academia Sinica. 

2 The source names W51A (G49. 5-0.4) or W51 are used synonymously in the following sections. 
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the envelope at a scale of 0.5 pc, resolving the eastern core in the JCMT map into the sources W51 
e2 and e8. W51 e2 is one of the strongest mm/submm continuum sources in the W51A region. In 
the highest angular resolution map obtained with the SMA at 870 um with 6 pa O '.'l , the polarization 
patterns appear to be pinched in e2 and also possibly in e8 hang et alJb009bh . In particular, the 
structure detected in the collapsing core e2 reveals hourglass-like features^ • These data currently 
provide the highest angular resolution information on the morphology of the magnetic field in the 
plane of sky obtained with the emitted polarized light in the W51A star formation site. A statisti- 
cal analysis based on a polarization structure function (of second order) shows a turbulent to mean 
magnetic field ratio which decre ases from the larger (BIMA) to the smaller (SMA) scales from 
~ 1 .2 to ~ 0.7 (|Koch et al.ll2010|) . possibly demonstrating that the role of magnetic field and turbu- 
lence evolves with sc a le. Collapse s i gnatu r es in W5 1 have been reported for various molecules in 



Rudolph et all dl 9901) :lHo & Yound(ll996l) : IZhang & Hoi d 19971) : I Young et alJdl998h : IZhang et al 
dl998h : ISolins et all toO^ . 



Figure [T reproduces the dust continuum Stokes / maps from the JCMT, BIMA and the SMA 
observations^, with their original (phase) centers at Right Ascension (J2000) a= 19 h 23 m 39 s .O, Dec- 
lination (J2000) 5=14°31W00; a= 19 h 23 m 44 s .2, 5=14°30'33."4 and a= 19 h 23 m 43 s .95, 5=14°30'34f'00. 
Overlaid in red are the magnetic field segments, rotated by 90° with respect to the originally de- 
tected polarization orientations. Typical measurement uncertainties of individual position angles 
are in the range of 5° to 10° and 3° to 9° for the SMA and for BIMA, respectively. The median 
uncertainty in the JCMT data is about 5°, with a few outliers between 30° and 40°. Additionally 
shown are the intensity gradient orientations in blue, which are relevant for the further analysis in 
the following sections. 

We remark that W51A was also observed with t he 350Atm polarime ter Hertz at the Caltech 
Submillimeter Observatory with a resolution of ~ 20" (IDotson et al.ll2010h . For our purposes here, 
we have found these data to be in agreement with and equivalent to the JCMT observation. They 
are, thus, not further discussed here. 



3. A Model Independent Approach: the Polarization - Intensity Gradient Method 

Various forces are interacting in molecular clouds. Maps of observed dust emission are re- 
flecting the overall result of gravity, pressure and magnetic field forces, and possible additional 



3 Regardless of the different physical scales in the JCMT, BIMA and the SMA observations, areas in a any map 
identified by clear emission peaks are denoted as cores, here and in the following sections. 

4 The reduced JCMT data are available at http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/ 
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constituents. Dust emission and magnetic field morphologies are left with a geometrical imprint 
by the combined effect of all these forces. Interestingly, polarization orientations are often ob- 
served to be close to tangential to dust emission intensity contours. Therefore, magnetic field and 
intensity gradient directions show a correlation, wh ere the difference 5 in their orientations can 



be linked to the magnetic field strength (Figure 3 in iKoch et alj (|2012l) ). In this method, various 
components in the ideal magneto-hydrodynamic (MHD) force equation are identified in dust po- 
larization and Stokes / maps. In particular, a change in emission intensity (gradient) is assumed 
to be the result of the transport of matter driven by the combination of all the forces in the MHD 
equation. Adopting this, it then follows that the gradient in the dust emission Stokes / intensity 
defines the resulting direction of motion in the MHD force equation. A force triangle, where the 
vector sum of all forces is set equal to the intensity gradient, can then be constructed (Figure 3 in 



Koch et al.l (120 12|) ). As a result, the total magnetic field strength as a function of position in a map 



can be calculated as: 



sin a 



where the angle ip is the difference in orientations between the gravitational pull and the intensity 
gradient, and the angle a is the difference between the polarization and the intensity gradient 
orientations. pV0 and VP are the gravitational pull and the pressure gradient, respectively. R is 
the magnetic field radius. Generally, all variables are functions of positions in a map. In the case 
of W51 e2, when neglecting the pressure gradient, the field strengths vary between ~ 2 mG and 



~ 20 mG with a radial profile B(r) ~ r 1//2 (|Koch et al.ll2012|) . The field strength averaged over the 
e2 core is ~7.7 mG. 

In order to determine the field strength B in Equation (QQ), the mass (gravitational potential <p) 
and the pressure gradient, if significant, need to be known. Contrary to this absolute field strength 
measure, the relative importance of the field compared to the other forces is directly imprinted 
in the field and intensity morphologies. With the magnetic field tension force F B = -f-^ and the 
gravitational and pressure forces \F G + F P \ = |pV0 + VP|, Equation (Q~|) is rewritten as: 



Sin ^' (2) 



\ F G+Fp\J i oca i Vsina //0ffl/ 

where we have introduced E B to define the field significance. The polarization - intensity gradient 
method, thus, provides a way to estimate the local magnetic field significance relative to other 
forces in a model-independent way. It is based on a generally valid (ideal) MHD equation, but it is 
independent of any molecular cloud/core models. Furthermore, this technique to extract the field 
significance is free of any necessity of mass and field strength measurements. The ratio in Equation 
© is solely based on measured angles which reflect the geometrical imprint of the various forces. 
Consequently, this also means that molecular clouds with accordingly scaled masses and field 
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strengths can show identical morphologies. This self-similar (or scale-free) property is lifted when 
calculating the field strength B for a particular cloud mass in Equation CD- 

In any case, the angle factor = S B provides a quantitative criterion as to whether the 
magnetic field can prevent an area in a molecular cloud fr om gravitational collapse (E B > 1) or not 
(S B < 1). As further demonstrated in iKoch et al.1 (|2012l) . T, B is only minimally or even not at all 
affected by projection effects. 



4. Results 

4.1. Magnetic Field - Intensity Gradient Correlation 

The dust continuum Stokes / maps of the JCMT, BIMA and the SMA are shown in the left 
column panels in Figure \\\ In this sequence of subsequently higher resolution maps, the JCMT 
main core (~ 40" in size) is resolved into two cores (~ 8") in the BIMA observation, where the 
Northern core is further resolved into W51 e2 (~ 2") in the SMA map. Overlaid on the dust 
continuum maps are the magnetic field segments (red) and the intens ity gradient segme nts (blue). 



For completeness, the SMA result for W51 e2 is reproduced from IKoch et al.1 (I2012|) . A tight 



correlation between magnetic field p osition angle (P. A. ) and the corresponding intensity gradient 



PA. (angle 5 < tc/2 in Figure 3 in IKoch et al.1 (|2012l) ) is obvious for many of the pairs. The 



correlation coefficients - in the definition of Pearson's linear correlation coefficient - are 0.71, 
0.72 and 0.95 for the JCMT, BIMA and the SMA data, respectively. Thus, there is seemingly a 
trend for an increasing alignment between magnetic field and intensity gradient orientations with 
smaller scales, with the tightest correlation found for the collapsing core. This correlation is being 
investigated in a separate work on a larger data sample (Koch et al. 2012; in preparation). It 
suffices to mention here that the above correlation can possibly serve as an indicator for the role of 
the magnetic field through the evolutionary stages of a molecular cloud. 



4.2. Relative Magnetic Field Significance from Large to Small Scales 

Maps of the magnetic field to gravitational force ratio, E B = 22^, are displayed in the right 
column panels in Figure [TJ The pressure gradient force, VP, typically being small compared to 
gravity, is omitted here. For the JCMT and BIMA data, two gravity centers at the two emission 
peaks in each map are assumed in order to calculate ip. The angle ip is then simply measured in 
between the intensity gradient direction (left column panels in Figure [U and the gravity center 
direction. For W51 e2 (SMA data), a single gravity center at the emission peak in the SMA map is 
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adopted. The angle a = n/2 — 5, where 5 is the difference in between the intensity gradient and the 
magnetic field orientations, is deduced from the left column panels in Figure [IJ All three maps - 
originating from three different instruments, one single dish and two interferometers - show clear 
differences in the ratios between core regions and areas in between the cores or further away. With 
the largest mapping area, the JCMT observation reveals very distinct features. With the exception 
of only two segments in the eastern core, both core regions show ratios below 0.5. In between 
the two cores the ratios reveal peaks in between 1.5 and ~5, with large stretches in the northeast- 
southwest direction larger than one. The BIMA observation is rather restricted to the cores, with 
an area in the northern core without polarization detection. Nevertheless, both cores clearly show 
ratios below one, except one segment in the southern core. The four segments in east- west direction 
(around y-offset ~ -4) reveal ratios larger or around one. Though limited to a few segments only, 
this still points toward a trend of increased ratios outside the immediate core areas. The detection 
by the SMA mostly resolves polarization features in the e2 core. The ratio averaged over the core 
is ~ 0.2. The north-west extension shows two segments with a ratio larger than one, and an average 
of about 0.8. If assuming an additional new core being formed here, the ratios get reduced to about 
0.5. In any case, the main collapsing core very clearly shows ratios below one, whereas the more 
distant north -west extensio n shows larger values. We remark that the total outflow mass in e2 



(ITang et al 



Shi et al.1 (|2010h ) is negligibly small compared to the total core mass of about 220 M & 
2009bh . Additionally, with a dust-to-gas ratio of about 1 to 100, the outflow is unlikely 



to be detected in the dust continuum with the SMA sensitivity. Consequently, we do not expect 
the magnetic field and the dust continuum Stokes / morphologies to be significantly affected by 
outflows. Our analysis, therefore, still quantifies the field to gravitational force ratio. Th e SMA e8 



core is not further analyzed here because only a few polarization segments are detected (|Tang et al 



2009b). 



In summary, the cores in the JCMT, BIMA and SMA data have average ratios of 0.33, 0.45 
(east, west), 0.38, 0.49 (north, south) and ~ 0.2 (main core), respectively. The other areas reveal 
significantly and systematically larger ratios ( £ 1). Despite probing three very different physical 
scales, the three data sets equally reveal a minor role of the magnetic field, F B < F G , in the center 
regions, and an increasingly more significant field, F B £ F G , at larger distances. 

The distinct features of the field-to-gravity force ratio in the right panels in Figure \T\ are 
further analyzed in the following. Azimuthally averaged radial profiles, binned at half of the beam 
resolution, are displayed in Figure[2] The profiles are centered at the dust continuum peaks. From 
the JCMT data, the main (eastern) peak covering the BIMA observation and the SMA e2 core is 
selected. On the next smaller scale, the BIMA southern core is chosen because its force ratio in 
the center has a smaller uncertainty than the northern core. As already manifest in the right panels 
in Figure [H a clear difference between center and outer regions is seen. Both the JCMT and SMA 
profiles show ratios across the center (~ 0.3 for JCMT and ~ 0.2 for SMA) with relatively little 
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variations. This is followed by a rather abrupt increase over a short distance (one bin) by a factor of 
~ 2 (JCMT) and ~ 4 (SMA). A plateau, with ratios ~ 0.8 for both JCMT and SMA, then extends 
over distances comparable to or larger than the center regions. The JCMT data, with the largest 
mapped area, then show another abrupt change to a ratio beyond one. No polarized emission is 
detected for W51 e2 at larger distances. The case of BIMA is less clear, but still shows values in 
the center of ~ 0. 1 and ~ 0.6, followed by a plateau around 0.6 and an increase to larger than one 
at the largest distances. For comparison, azimuthally averaged emission intensity profiles are also 
shown. 

For completeness, profiles for the JCMT western and the BIMA northern core are shown 
in Figure |3] The result for the very central region in the BIMA northern core is less conclusive 
because only a single value around 0.7 with a ±50% error is available for the force ratio. Never- 
theless, values are clearly below one in the core region with a trend to grow at larger distances up 
to about 0.9. Similarly to the eastern core, the western core in the JCMT data shows ratios around 
0.2 to 0.6 up to about 20" in distance. Further away from the peak, ratios grow up to about 5 with 
some oscillatory behavior remaining larger than one. Thus, despite being less pronounced, the two 
additional cores here also reveal a transition in their force ratios between central and more distant 
ratios. In particular, in the overlapping region between the BIMA northern core and the SMA e2 
core, the ratios show similar trends regardless of the different field morphologies (hourglass-like 
for SMA, more uniform for BIMA). 

Errors in the field-to-gravity force ratio are calculated by propagating the measurement un- 
certainties Atp and Aa through Equation ©. The measured magnetic field PA. uncertainty, in 
the range of a few degrees to ~ 10°, is assumed for Aa. A typical uncertainty of 3° -5° in Aip 
results after interpolation when calculating the intensity gradients. This leads to errors for individ- 
ual ratios in Figure[2]of ~ ±4% to ±19%, ~ ±4% to ±40% and ~ ±2% to ±20% (with a single 
outlier at ±60%) for the SMA, BIMA and the JCMT data, respectively. After averaging in each 
bin, the errors are typically reduced due to the sample variance factor, resulting in average errors 
of ~ ±10% or less. Average errors of ~ ±20% remain for two bins in the JCMT data. A 10% 
measurement uncertainty is conservatively estimated for the emission intensity when calculating 
errors for its radial profiles. Binned errors are then typically at the percent level (Figure [2]). Except 
for the single central value for the force ratio in the BIMA northern core (~ ±50% error), similar 
errors are present in Figure |3j 

Finally, it is important to keep in mind that the correlation presented in Section [47T1 (based on 
one angle between two orientations) can be affected by projection effects. Generally, all values 
are integrated along the line of sight. The field significance S B presented he re is much less o r not 



at all affected by projection, because it is based on the ratio of two angles (|Koch et al]|2012|) . It, 
nevertheless, still deals with quantities averaged along the line of sight. 
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5. Discussion 

5.1. Self-Similar Profiles 

Given the common features in the radial profiles of the magnetic to gravitational force ratio 
(Figures [2] and [3]), we address here the question of self- similarity. We focus on the JCMT main 
core where the BIMA and SMA data provide higher resolution follow-up observations. The BIMA 
southern core is adopted due to its better statistics. All three data sets in Figure |2] show a rather 
constant ratio across their cores (plateau-like or a single binned radius), followed by a plateau 
with a larger ratio before eventually the ratio grows to values larger than or around one. With 
the cores in the maps of Figure Q] being clearly identified, we choose to normalize the distances 
from the peaks to units of core sizes; i.e. all profiles are aligned to one normalized core size. The 
normalization radii, chosen to be those bins where the emission intensity profiles in Figure [2] flatten 
out, are 21."85, 3''42 and l."06 for the JCMT, the BIMA and the SMA data, respectively. At these 
radii, the emissions are ~ 15% to 20% of the peak emissions. Figure @] shows the result. 

Even with some uncertainty left in the normalization, it is obvious that the aligned profiles 
show a close resemblance. This self- similarity suggests that the interplay between magnetic and 
gravitational force is independent of the three different scales probed here. Thus, the analysis here 
quantifies the relative magnetic field significance from being generally dominant or comparable to 
gravity at distances twice the core size to minor or negligible inside the core. We stress that this 
result seems to hold generally for all the different field morphologies analyzed here, including even 
the more irregular cases shown in Figure [3j 



5.2. Mass-to-Flux Ratio Derived with Polarization Intensity- Gradient Method 



In this section we investigate the connection between the force ratio in Equation Q and the 
mass-to-flux ratio for a molecular cloud. The mass scale ass ociated with the amount of magnetic 
flux $ threaded by a self-gravitating cloud was introduced in LMestel & Spitzeil (| 19561) . The mag- 



netic critical mass M$ = 2nG i /2 , with $ = nR 2 B where R is the cloud radius, defines the maximum 
mass that can be supported by th e magnetic field against gravitational collapse if no other forces 
are present (e.g. IShu et al.1 (|1999m . This leads to the notions of magnetically supercritical clouds 
when M > M$, and magnetically subcritical clouds when M < M$. Contrary to the mass-to-flux 
ratio - which is typically applied to the entire (global) molecular cloud - the force ratio in Equa- 
tion © is a local criterion comparing the field significance with gravity at a specific location in 
a cloud. Neglecting pressure gradients and explicitly writing out the local field force and gravity 
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force terms, it can be expressed as a ratio of local flux over local mass: 



smipe _ $g nu J_ J_ r2 (3) 

sina £ ~ mj ' M t ' R B ,e Rj ' G/ ' 4vr 3 G' 
where the lower index t refers to a local quantity. Mi is the gravitating mass leading to the local 
gravitational force F Gi = acting upon a local mass element m e at a distance R G t . is the 

K G,l 

flux associated with the local mass within a flux tube of radius Ri (Figure [5]). Rb c is the local field 
radius. G and rj are the gravitational constant and a numerical unit conversion factor, respectively. 
Equation © is generally valid. However, in order to make further use of it, the local mass element 
mi and Mi need to be specified. We, therefore, make the basic assumption of spherical symmetry. 
The enclosed mass within a distance r from the center is then Mi(r) = ^Anr' 1 p(r')dr' = N r rhi{r), 
where N r is the number of local flux tubes with an average local mass rhi(r) where both depend on 
the integration radius r. A local mass-to-flux ratio then results from Equation ©: 

«« _ / sin* y* / •» f j^yn 

(pi \smaij \mt(r)J * \GJ 

R 2 

where we have used M e (r) = N r m~i(r) with Af r = -Jr. We note that the above expression is still fairly 
general and valid for any density profile p(f). In particular, since the ratio depends on the ratio of 
the mass profile over a mean mass, only the functional form of the density profile is relevant, but not 
its absolute value. Introducing a mean local mass, m~e(r), allows us to write the local mass-to-flux 
ratio with the explicit spherical radius dependence only in the second term on the right hand side 
of Equation (HI). Besides the spherical symmetry assumption for the mass distribution, Equation 
© is basically valid at any position in a molecular cloudj. The average inverse of the force ratio 

/ \ 1/2 

decreases with larger radius (FigureH]), and ( ) < 1 is monotonically decreasing for centrally 



peaked density profiles, with unity in the center. The local mass-to-flux ratio in Equation © 
is, thus, decreasing with larger radius from the center. Azimuthally averaging Equation © and 
normalizing it to the critical mass-to-flux ratio, we get the local normalized mass-to-flux ratio - for 
individual subvolumes of fixed size R e in a cloud - as a function of radius: 

/AM \ //sin^V 1/2 1/2 \ (mi{r)\ 12 R? _ 1/2 / /sin^A " 1/2 \ _ 1/2 

(5) 



5 We note that it is also possible to proceed directly with Equation (0. Assuming that both mi and Mi are pro- 
portional to the integrated dust emission with an identical conversion factor, the ratio ^ can directly be calculated 
from a dust continu um emission map. In this way, it is valid for any arbitrary cloud shape. This is identical to the 
approach outlined in lKoch et al.l (120121) where a local gravity direction is derived for any cloud shape. Equation (0), 
in its most general form, then leads to a map of local mass-to-flux ratios. However, we are here aiming at revealing 
changes in the mass-to-flux ratios with radius. Therefore, the spherical symmetry assumption and azimuthal averaging 
are appropriate (Equation (O). In practice, we will calculate ^ from the dust emission maps as mentioned above. 
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where < ... > r denotes azimuthal averaging at ra dius r. When norm alizing to a local flux tube, we 



can simply set R = R e . Furthermore, as found in lKoch et al.1 (I2012|) . R BI is roughly constant over 



an observed map, and similar to Ri, which we identify with the beam resolution of an observation. 

R ll2 R- { 

With this simplification, BJ _ 3 / 2 e ~ 1 • Assuming the cloud gravitating mass to be proportional to the 

R 

dust emission, ^ is derived from the observed emission intensity profiles in Figure [21 As already 
indicated with the notation, Equation © quantifies the differential mass-to-flux ratio as a function 
of radius. 

In a next step, we derive an integrated mass-to-flux ratio, i.e. we are asking whether the entire 
cloud inside a certain radius r is magnetically supercritical. We then have to evaluate: 

M T Nr m e 
<<r)=^ r -^=M l (r)- > > • (6) 




(fie can be expressed with Equation ©. The summation is over all local mass and flux elements 
within the radius r, where the summation limit N r depends on r. After some algebra and after 
normalizing with the critical mass-to-flux ratio we find: 

$ J norm \ \^Ot t J ' / ^ m(j) J r 1 1 2 

where the averaging < ... > n is for each bin r,, with N n being the number of subvolumes within 
ri and r i+ \. Nb is the number of bins within r (Figure [5]). Equation © states that the integrated 
mass-to-flux ratio is calculated by adding bin-averaged force ratios which are weighted with a 
mass function and N r ., where the latter one results from the integrated volume growing with r. 
The normalization is with respect to r, i.e. the increasingly larger cloud (with growing radius r) is 
normalized with its corresponding critical mass-to-flux ratio. Ri is constant and again set by the 
beam resolution. R B £ is also roughly constant. When omitting the summation, setting N n = 1 and 
r = R e and replacing m(f) with m(r,), Equation © reduces to the differential mass-to-flux ratio in 
Equation ©. 

Finally, a global mass-to-flux ratio for the entire cloud as a single entity can be derived from 
Equation © in the limiting case where all (small) local quantities grow to cloud-size: mi — > M 
and < R B e >~< Re >~< Rg >= R- The force ratio needs to be averaged over the entire cloud: 
— > ( 2HU££ \ This then leads to the global mass-to-flux ratio, normalized to the critical mass- 



sincy \ sina£ 

to-flux ratio with R 



M\ / / sin^A 1 \ ( R 



$ / \ \ sina* / / \Ro 



!/2\ / n \ "3/2 



vr- 1 / 2 , (8) 
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where Rq = R is the cloud radius. 



The set of Equations ©, and ([8]) defines a differential, integrated and global mass-to-flux 
ratio, respectively, which can readily be calculated from observations. These equations provide 
the criteria to evaluate whether a subvolume (local or integrated up to a certain radius) or an entire 
cloud is magnetically supercritical. In all cases, the inverse square root of the force ratio is the 
key measure leading to an estimate independent of any mass or field strength input. Except for 
the global mass-to-flux ratio, the local ratios need the functional form of the mass profile as an 
additional input. 

A schematic illustration of the different subvolumes analyzed by the various ratios is shown 
in Figure[5J Figure[6]illustrates these mass-to-flux ratios for the SMA observation. The differential 
mass-to-flux ratio shows a clear transition from a magnetically subcritical state at larger radii to a 
supercritical state at smaller radii. This might indicate that a fixed subvolume can only collapse 
closer to the center, but not in the outer parts of a cloud/core. This finding possibly also provides 
an explanation for why fragmentation and clustering of stars preferentially occur at the center bu t 
not at larger radii. Such a result is further consistent with the theoretical work by lShu & Lil(|1997l) . 
where the dilution of self-gravitational forces by the magnetic tension force (see Section 15.31 and 
Figure |U) is found to reduce the tendency of a disk to fragment. As further utilized in Section [531 
this dilution is more significant at larger radii. Consequently, the fragmentation process is more 
significantly suppressed here. The number of individual collapsing subvolumes in the center region 
(with smaller or non-existent dilution) might then be linked to the multiplicity of s tars. Addition- 
ally, the finding here might give support to the suggestion by iMouschoviasI (I1976Q that magnetic 
tension might be able to hold back th e envelope and, t herefore, prevent it from following the cloud 
core collapse. In a related later work. IShu et al.l (|2004|) investigated whether magnetic tension can 
define the masses of forming stars by holding up the subcritical envelope of a molecular cloud 
while its supercritical core is collapsing. A split monopole - as it is possibly seen in the field 
morphology in the e2 core (ITang et alj|2009bl) - leads to magnetic levitation in their work. How- 
ever, they find that magnetic suspension alone can not keep the subcritical envelope from falling 
onto the center, but additional dynamic levitation by winds and jets is needed for a halt of infall. 
The results here are also in qualitative agree ment with the recent numerica l simulation of giant 
molecular clouds with ambipolar diffusion by Vazquez- Semadeni et al. ( 201 1|) . They find that the 
mass-to-flux ratio inside a cloud exhibits large spatial fluctuations, spanning an order of magnitude 
or more. Thus, local clumps can become magnetically supercritical within a globally subcritical 
system. Esti mating 1Z (the ratio be tween the mass-to-flux ratio in the core and in the envelope as 
introduced in lCrutcher et al.l (|2009|) ) from the values at the largest distances and the values close to 
the center (Figure [6]), we find 1Z > 1 . This would be in favor of ambipolar diffusion driven mod- 
els. However, the recent statistical a nalysis of 1Z based o n numerical simulations with supersonic 
magnetohydrodynamical turbulence (Bertram et al.ll201 if) yields a large scatter in 1Z with values 
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both smaller and larger then unity. It, thus, has to be acknowledged that the interpretation of 1Z 
as a measure to distinguish between ambipolar diffusion and turbulence driven star formation is 
possibly not unambiguous. 

The integrated mass-to-flux ratio in Figure |6]is driven by the increasing volume. It, therefore, 
grows with radius. As expected, this ratio is generally larger than the differential ratio of a subvol- 
ume, because mass scales with volume whereas flux scales with area. The slight mismatch in the 
innermost bin results from different weightings for the two ratios. Finally, the global mass-to-flux 
ratio is close to the integrated ratio at large radii. The difference here results from the incomplete 
sampling of N n in azimuth for those bins which contain some depolarization zones. It is important 
to remark that both the integrated and global ratios present average values for the entire cloud. 
Therefore, they are likely biased toward some subvolumes. Whereas these ratios provide a de- 
scription as to whether a cloud as an entity is in a sub- or supercritical state, they can not properly 
assess the state of individual subvolumes. This question can only be further addressed with the 
differential mass-to-flux ratio. The finding here clearly demonstrates the differences and the need 
for a measurement of the local field significance. 



5.3. Dynamical Implication and Star Formation Efficiency 

In this section we aim at investigating further implications of the local force ratio |^ (Figure 
Bi and the differential mass-to-flux ratio (Figure©. We limit the discussion here to the collapsing 
core W5 1 e2. Global and/or average cloud properties are typically used to characterize its state; e.g. 
in order to determine whether a system is gravitationally bound, or whether it can collapse or not. 
With the local criteria - the magnetic field to gravity force ratio in Equation © and the differential 
mass-to-flux ratio in Equation © - we have tools to go one step further. We can not only address 
the question whether a cloud collapses or not, but we can even ask: Where does a cloud collapse? 
How does it collapse? Does all the gas take part in the collapse? And finally, depending on the 
answers to these questions, what are the consequences for the inferred star formation efficiency? 

The change with radius in the magnetic field to gravity force ratio (Figure clearly demon- 
strates that gravity is not everywhere equally efficient to initiate and to keep driving a collapse. 
The result in Figure @] suggests and quantifies an effective gravitational force which is reduced by 
the presence of the magnetic field. I ndeed, on t he the ory side, a concept of diluted gravity (by the 



magnetic field) was put forward by IShu & Lil (|1997D . Thus, we can define an effectively acting 
gravitational force as 

F*(r) = f 1 - ) . Fmax(r) =< ec >r .F max (r\ (9) 



sin a 

where < ... > r denotes azimuthal averaging at radius r and F max (r) is the maximum non-diluted 
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gravitational force resulting from an enclosed spherically symmetric mass distribution within r. 
is only defined where < 1 , which possibly holds only for a limited area. In Equation © we 
have introduced the gravity efficiency factor e G G [0,lfl. " 

Besides the diluted gravity which slows down the accretion and collapse process, the fraction 
of volume of a cloud (or core) that actually can collapse is relevant; i.e. where is the cloud/core 
magnetically supercritical? The differential mass-to-flux ratio (Figure© provides a local criterion 
that answers this question. The increasing ratio toward the center indicates that only the central part 
of the core has accumulated enough mass to overcome the magnetic flux. By reading off the radius 
where the core state changes from magnetically sub- to supercritical, the potentially collapsing 
volume and its associated mass can be estimated. This defines a volume efficiency, e v G [0, 1], as 
compared to the maximum efficiency e v = 1 where the entire core can collapse. 

Based on the two above described findings, we proceed to estimate a star formation efficiency 
with the assumptions: (1) the available effective gravitational force is diluted by the presence of 
the magnetic field and (2) only the regions where the differential mass-to-flux ratio is larger than 
one will eventually collapse and form stars. We choose to reference our estimate to the well-known 
pressure-less free-fall collapse of a spherical gas sphere. The dynamics in this case are governed 
by the momentum equation ^ = — G^f, where the convective derivative is ^ = ^ + for the gas 
infall velocity v. The enclosed mass within radius r is M r = j^Aur' 2 p(r')dr' with the gas density p. 
With the diluted gravity, the momentum equation reads: 



Dv GM r 
— =-<e G > r — 
Dt r 



<e G > r ~. (10) 



Technically, the gravity efficiency factor < e G >,- could be absorbed into an effective density p*(r). 
Therefore, the general solution (for any M r ) of the Equation (fTOl) for a fluid element starting at a 
distance r and reaching the center at r = at the time t is still valid. The time U in the presence of 
the magnetic field is theifl: 



TV 



U = - x \ 1° ^<e G >- l ' 2 -t (11) 

2]j2GMMr),<eG>r) 



6 Unlike the gravity efficiency factor ec, the field significance Eg in Equation (f2]) is not constrained to an upper 
limit of one. 

7 For a uniform density and < tQ > r = 1, Equation (fTTI) describes a synchronized collapse with the well-known 
free-fall time ?// = y^ff^;. More realistic density profiles with a central concentration, e.g. p(r) = p c (jr^J lead to 

r— / \ a/2 

a collapse time t = J ^""^ ■ iy-J which increases with larger distance ro from the center and then reproduces an 
inside-out collapse. 
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Strictly speaking, < e G >r in the above equation is still a function of radius. In writing the right 
hand side we have assumed it to be averaged within an appropriate radius. The average accretion 
rate < M >, assuming the entire mass M within a cloud/core to be accreted and to collapse, is 
< M >= M/t. Adding the volume efficiency factor e v together with Equation CCD, the effective 
accretion rate < M* >, corrected for our observed magnetic field features, becomes: 

»v M * e vM 

U < ec >" 1/2 -t 

The resulting star formation efficiency relative to the standard free-fall one is then expressed as: 

-t— «ey<e G > / • (13) 

<M> 

We note that the above phenomenological modeling is still valid in the presence of additional force 
terms in the momentum equation with any density profiles. The pressure-less free-fall collapse 
simply allows us to analytically express the results. Figure [7] shows the reduced star formation 
efficiency taking into account the effects of the magnetic field. The red dots indicate efficiencies 
estimated from the gravity dilution and the reduced collapse volume based on the Figures |4] and 
|6j Since 1 - ^j"^ ~ at large distances, an arbitrarily low star formation efficiency would result 
from adopting these values for e G . Instead, we adopt values averaged over the entire core (force 
ratio ~ 0.33), the central core only (~ 0.2) and the outer plateau (~ 0.85). We consider these 
values adequate to estimate and sample the gravity dilution for W51 e2. Similarly, for the volume 
efficiency factor we read a radius of about 1" within which the core is magnetically supercritical 
as compared to an overall size of about 2". For comparison, a larger volume of about 1.5" is 
considered as well. The resulting efficiencies are reduced at least to a conservative value of ~ 0.35 
or less when still assuming the larger volume to collapse. The efficiencies drop to ~ 0. 1 with 
the smaller volume, with a likely limit of about ~ 0.05 in combination with the largest observed 
gravi ty dilution. These values are close to the observationally inferred efficiencies of a few percent 

I e.g., Krumholz & Tanll2007|) . For comparison, recent numerical simulations by iNakamura & Li 
201 1|) find a star formation efficiency per global free-fall time of ~ 20 to ~ 30% depending on 
the magnetic field strength. When they further take into account the feedback from outflows, their 
efficiencies can be reduced by another order of magnitude. 

We finally remark that the discussion here is limited to the core e2. A more complete picture, 
estimating the star formation efficiency starting from the largest scales of the cloud envelop, will 
need to dynamically link the observed self-similar profiles from Figure 51 Caution is needed here 
because linking structures of different size scales is non-trivial due to additional structures possibly 
generated by, e.g. fragmentation and turbulence. Nevertheless, this first estimate further solidifies 
the significance and impact of the magnetic field. Moreover, it demonstrates that the magnetic 
field is capable of reducing the star formation efficiency by at least one order of magnitude or even 
more. 
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6. Summary and Conclusion 



We have applied the polarization - intensity gradient method (IKoch et al.H2012|) to a set of 
single dish and interferometer dust continuum data of the W5 1 star formation region. These obser- 
vations cover scales from the initial parent cloud (~ 2 pc, JCMT) to the large-scale envelope (~ 
0.3 pc, BIMA) and to the collapsing core (~ 60 mpc, SMA). Besides leading to a magnetic field 
strength as a function of position in a map, our new method also provides a way to estimate the 
local magnetic field to gravity force ratio. The technique is model-independent, without any input 
from mass or field strength, solely making use of measured angles in dust polarization and Stokes 
/ map. Here, we have focused on this force ratio and on some of its implications. Our main results 
are summarized in the following. 



1. Magnetic field - intensity gradient correlation: The correlation in the magnetic field and 
dust intensity gradient orienta tions - which served as a starting point for the new method 
developed in lKoch et al. ( 2012 ) - is also found at larger scales in the BIMA and JCMT data. 



The correlation coefficients seem to be the larger the higher the resolution is; i.e. the tightest 
correlation is found for the collapsing core e2 (Figure [H left panels). 

2. Magnetic Field Significance: Maps of the magnetic field to gravity force ratio show distinct 
features for all three data sets (Figured! right panels). At larger distances from the emission 
peaks, the field force is comparable or larger than the gravity force (ratio ^ 1). Closer to 
the emission peaks, gravity becomes more dominant (ratio < 1). This suggests and quanti- 
fies that gravity is not everywhere acting equally efficiently, but is being more significantly 
opposed by the magnetic field tension force at larger distances. Based on this finding, we 
introduce a gravity efficiency factor. This also seems to be in agreement with the concept of 
gravity dilution (FigureH]). Moreover, this is in support of inside-out collapse scenarios. 

3. Self-similarity: Azimuthally averaged radial profiles of the force ratios show similar features 
for the JCMT, BIMA and the SMA data (Figure [2]). When normalized to core sizes, these 
profiles closely align (FigureH]). This points toward self-similar properties in the magnetic 
field and gravity interplay from the large parent cloud down to the collapsing core. 

4. Mass-to-flux ratio: The force ratio can be converted into a mass-to-flux ratio. In its most 
general form, the mass-to-flux ratio can then be expressed as the inverse square root of the 
force ratio multiplied by the square root of a ratio of two mass elements. With the force ratio 
being a local criterion, this naturally leads to a local (differential) mass-to-flux ratio for a 
subvolume in a molecular cloud as a function of position. An integrated and global mass- 
to-flux ratio can be derived by appropriately adding and averaging local quantities. Similar 
to the force ratio, all the various mass-to-flux ratios do not have to rely on absolute mass 
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and field strength inputs, but they are largely model-independent. This finding then does 
not only provide a criterion to decide whether a molecular cloud as an entity is magnetically 
supercritical or not, but it extends this criterion to any subvolume (Figures [5] and [6]). In the 
case of W51 e2, a transition occurs from subcritical at larger distances to supercritical in the 
central area. This might explain why fragmentation is found to happen preferentially in the 
center. 

5. Star formation efficiency: A reduced star formation efficiency (compared to the non-magnetized 
free-fall case) is derived based on two observed magnetic field properties: (1) Gravity dilu- 
tion (Figure H]) defines an effectively available gravity force to initiate and drive a collapse 
and, therefore, it increases the collapse time. (2) The differential mass-to-flux ratio (as a 
function of radius) provides a hint that the magnetic field reduces the volume of supercriti- 
cal state and, therefore, only a limited volume can collapse and form stars (Figure©. With 
these two findings a modified accretion rate can be calculated. Comparing this to a standard 
free-fall rate shows that the magnetic field is capable of reducing a standard star formation 
efficiency to 10% or less. 

The authors thank the referee for valuable comments which led to further important insight in 
this work. P.T.P.H. is supported by NSC grant NSC97-21 12-M-001-007-MY3. 
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Fig. 1 . Left column panels: Dust continuum Stokes / o bservations towa rd W51. From top to bottom: JCMT a t 850 jan (9 ft 9" 3 , from 

IChrvsostomou et alj 120021) '). BIMA at 1.3 mm (0 a 3"0, from lLai et all l200ll) ) and SMA at 870 fim (d ft 0"7, from lTang eta?] J2009bl) ). The 
beam resolutions are shown with the black ellipses in the lower left or right corners. The BIMA observation is zooming in onto the main core 
(dashed rectangle) in the JCMT map. The BIMA main core (dashed square) is further resolved in the SMA map. The cores e2 and e8 are labeled. 
Colors correspond to the color wedges on the right hand side with units in Jy/beam. Overlaid are the magnetic field segments (red segments) and 
the intensity gradient segments (blue segments) at the locations where polarized emission was detected. Magnetic field segments are plotted by 
rotating the polarization segments by 90°. The difference in P.A.s between the magnetic field and the intensity gradient is the angle 8 (Section 
I4.lt . The length of the segments is arbitrary, and for visual guidance only. The axes offset positions (in arcsec) are with respect to the original 
(phase) centers of the observations. At the distance of W51, 1" is equivalent to about 30 mpc. Right column panels: The ratio of the magnetic field 
force Fb compared to the gravitational force Fo, £b = jjjj^j, as a function of position. The pixel scale is identical to the gridding in the left column 
panels. From top to bottom are shown the maps for the JCMT, BIMA and the SMA observations. Overlaid in black are the Stokes / dust continuum 
emission contours from Figure [T] Colors correspond to the color wedges on the right hand side, displaying the dimensionless ratio. Blue colors 
indicate where gravity dominates over the magnetic field (Fg < Fq). 
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Fig. 2. The magnetic field to gravitational force ratio, Sb = as a function of distance from the emission peaks (blue). From top to 

bottom shown are the JCMT, BIMA and the SMA results. The corresponding intensity emission profiles are displayed in green with the axis at the 
right hand side. Values are azimuthally averaged and binned at half of the beam resolution. From the JCMT and BIMA observations, the results 
for the main peak and the southern peak, respectively, are shown. Binned values are connected with lines for visual guidance only. Errors for the 
force ratios are calculated by propagating typical measurement uncertainties in rjj and a through Equation |2). A 10% uncertainty in the intensity 
emission is assumed. Error bars at the smallest and largest scales are typically growing because of a smaller sample variance factor. Binned errors 
are all around ~ ±10% or less, with the only exception being the JCMT data at the largest scale with an error of ~ ±20%. 



-22- 



W51 - JCMT 850 p.m: western core 

6 1 1 1 1 1 1 1 60 




distance from peak [arcsec] 
W51 - BIMA 1 .3 mm: northern core 



g 

CO 



■> 0.5 
CO 
O) 

I 

o 

*T 

CD 



0.5 1 1.5 2 2.5 3 3.5 4 4.5 5" 

distance from peak [arcsec] 

Fig. 3. — As in Figure [2l the magnetic field to gravitational force ratio with the corresponding 
intensity emission profile as a function of distance from the emission peak for the JCMT western 
and the BIMA northern core. 
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Fig. 4. — The magnetic field to gravitational force ratio, as in Figure [2l but with the distances 
from the peaks normalized to units of the core sizes. With respect to Figure [2l the normalization 
radii are 21. "85, 3. "42 and I'.' 06 for the JCMT, the BIMA and the SMA data, respectively. The 
force balance = S B = 1) is indicated with the horizontal dashed black line. The effectively 
acting gravitational force (as introduced in Equation Q) is derived with the gravity efficiency 
< e G > r = 1 - (|^) r , which is symbolized through the lengths of the down-arrows. Similarly, the 
gravity dilution (by the magnetic field) is indicated with the dashed up-arrows. The error bars, 
identical to the ones in Figure [2l are omitted here. 
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Fig. 5. — Illustration of the various mass-to-flux ratios as derived in Section I5?2l Every small 
circle of radius Ri represents a local mass-to-flux ratio Equation ©). Its size is determined 
by the beam resolution. The differential mass-to-flux ratio (^(r), Equation ©) at radius r is 
calculated by azimuthally averaging all the local ratios of the hatched red circles inside the ring 
at radius r. The four solid red circles symbolize differential ratios at four different radii. The 
integrated mass-to-flux ratio (¥(< r), Equation ©) within radius r is shown with the blue hatched 
circle. Concentric blue circles with growing radii symbolize integrated ratios for increasingly 
larger volumes of the cloud. Azimuthally averaged ratios for each radius are displayed in the 
profiles in Figure [6J A single global ratio Equation ®) is calculated from the black circle 
encompassing the entire cloud. 
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Fig. 6. — The various mass-to-flux ratios normalized to the corresponding critical mass-to-flux ratios (Section l5^2l 
for W51 e2. The ratios are shown as a function of distance from the peak normalized to the core size. The solid 
line displays the normalized azimuthally averaged differential mass-to-flux ratios (Equation ((5), converted from the 

magnetic field to gravitational force ratios (Figure |2} with the additional mass ratio, ( j < 1. The mass ratio is 
derived from the integrated emission profiles in Figure[2] A clear transition between magnetically supercritical in the 
center to subcritical at larger distances is revealed. Open circles correspond to the individual ratios before averaging 
and before multiplying with the mass ratio. Thus, they are local upper limits (Equation (0). The dashed line is the 
integrated mass-to-flux ratio (Equation (O) normalized to the critical mass-to-flux ratio with its radius corresponding 
to each binning radius. The single global mass-to-flux ratio (Equation ([8])), averaging all force ratios and normalizing 
to the core radius, is displayed with the filled circle. Indicated with black dashed lines are the critical mass-to-flux ratio 

/ . \ -1/2 

and the core size (normalized to one). When propagating uncertainties through ~ ( ) , errors are suspect 

/ . \ -3/2 

to an additional factor \ I j as compared to the force ratios (Figure |2}. Thus, at larger radii where the force 
ratio is around one, errors in the mass-to-flux ratios are similar to or smaller than the ones for the field-to-gravity force, 
i.e. £ ± 10%. As the force ratios drop in the center area, errors grow larger by a factor of about 2 to 3, i.e. reaching 
~ ±20% to ~ ±30%. For clarity, error bars are not displayed here. 
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Fig. 7. — Estimated reduced star formation efficiency (contour levels) with a diluted gravitational 
force (efficiency factor eo) and a reduced collapsing volume (efficiency factor e v ) following Equa- 
tion (fl"3T) . The lower panel shows the zoomed-in area marked with dashed lines in the upper panel. 
The red dots are possible estimated efficiencies based on the force ratio (Figure H]) and the differ- 
ential mass-to-flux ratio (Figure [6]) for the case of W51 e2. 



